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Abstract 

We construct a model of wealth distribution, based on an interactive multiplica- 
tive stochastic process on static complex networks. Through numerical simulations 
we show that a decrease in the number of links discourages equality in wealth dis- 
tribution, while the rewiring of links in small- world networks encourages it. Inequal 
distributions obey log-normal distributions, which are produced by wealth cluster- 
ing. The rewiring of links breaks the wealth clustering and makes wealth obey the 
mean field type (power law) distributions. A mechanism that explains the appear- 
ance of log- normal distributions with a power law tail is proposed. 
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With the advance of the theory of complex networks [T|], it has become clear that many 
social networks can be understood as small-world networks @. On the other hand, recent 
high-precision studies of economic phenomena, helped by the availability of high-quality 
data in digital format ||, is starting to reveal the basic characterics of wealth distribution, 
which form a basis and at the same time a final product of economic activities. In view 
of these developments, we propose a stochastic model of wealth distribution built on a 
complex network and present the results of our numerical simulation in this paper. 

Several models for wealth distribution have been proposed with some reality and 
success [|], [5], [], [7|. They, however, lack one feature we consider important: In some 
models, each agent (economic body) goes through a stochastic process of increasing and 
decreasing wealth, quite independently from other agents. In some others, each agent 
interacts with either randomly-selected agents or neighboring agents on a fixed lattice. 
Apparently, such a selection of economic partnership/competition is far from reality, as 
is evident from complex network studies. This is one reason we propose our model, 
which features both the stochastic nature of economic activities and underlying network 
structures that are expected to be close to those of the relevant social activities. We 
also note that studies of economic phenomena that build upon the underlying complex 
networks are novel. 

In considering the mechanisms behind wealth distribution, it is important to include 
at least two basic natures of the wealth. One is a random change of the value of wealth 
and the other is the liquidity of it. As a minimal model that contains these features we 
chose an interactive multiplicative stochastic process || , which is defined by the following 
evolution equation: 

Wt(t + 1) = Oi(t)Wi(t) + a 2 W t {t) + J2 JiAtWAt) ~ E Jji(tWi(t), (1) 

where Wi(t) is the wealth of the i-ih agent at the time t, Oj(t) is a Gaussian random 
variable with mean m and variance 2a 2 , which describes the spontaneous increase or de- 
crease of the wealth. The rest of the terms in the r.h.s. of Eq. (H) describe transactions 
on business networks. We note that in Ref. || the process corresponding to Eq. (H) is 
described by a stochastic differential equation in the Stratonovich sense, which generates 
the second term proportional to a 2 when discretized. In addition, although naive dis- 
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cretization brings in a time interval At in Eq.(|T|), it has been absorbed by the rescaling 
of the parameters. 

The transaction matrix J^- is chosen based on the underlying network structure: If 
the i-th agent is directly connected to the j-th agent by a link, we take = J/Zi, where 
Zi is the number of agents connected to the i-th agent ("neighbors"). Otherwise, Jy = 0. 
The parameter J is a positive number, which represents the ratio of transactions and the 
wealth. Thus, the third term in the r.h.s. of Eq. (|I|) describes the incoming wealth, and 
the fourth term describes the outgoing wealth. Under this selection of transaction matrix, 
Eq. ([[]) can be written as follows: 

Wi(t + 1) = { ai (t) + a 2 - j} Wi(t) + JWi(t), (2) 

where Wi(t) is the mean wealth of the neighbors of the i-th agent. 

We consider one-dimensional and static network of iV-agents with periodic boundary 
conditions |], 0, which are determined by two parameters, i.e., a rewiring probability p and 
a mean number of links per agent Z. (We will also use a notation z = Z/N.) The rewiring 
probability p is defined by p = (the number of rewired links) / (the total number of links) e 
[0, 1]. At p = we have a regular network with Z{ = Z, in which the i-th agent is con- 
nected only with {i — Z/2, • • ■ , % — 1, i + 1, i + Z/2}-th agents. For p ^ 0, links are rewired 
to conserve the total number of links. At p = 1 we obtain completely random networks. 

In the mean field (MF) case, i.e., z ~ 1, the stationary probability density function 
p(w) is known to be the following in the large iV limit 0, 

n(m ) - (^-iyexp{-(/i-l)M 

where Wi(t) is the wealth normalized by the mean wealth W(t); Wi(t) = Wi(t) /W(i). (We 
will simply call W{(t) the normalized wealth hereafter.) The exponent is given by /i = 
1 + J /a 2 . In the large w range Eq. (|3|) becomes a power law distribution, p{w) oc w~^ 1+ ^\ 
and therefore the exponent \i is the Pareto index in economics. 

We have carried out numerical simulations for N = 10 4 , m = 1, J = 0.01 and 
cr 2 = 0.01. Under this parameter set, the Pareto index in the MF case is /x m f = 2. 
Some of the distributions we obtained are shown in Figs. [I] (a)-(d). In these figures the 
horizontal axis is the logarithm of w and the vertical axis is that of the rank. We have 
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Figure 1: Some of the wealth distributions obtained by numerical simulations. The rank 
distribution in regular networks is plotted in (a) for z — 0.1 and in (b) for z = 0.001. The 
rank distribution in small-world networks is plotted in (c) for z = 0.001 and p = 0.01 and 
in (d) for z = 0.001 and p = 0.5. 

numerically calculated up to 2 x 10 4 time steps, and drawn distributions at each 4 x 10 3 
step. The rank distributions in the case of regular networks with z = 0.1 correspond to 
Fig. |T] (a). The solid lines show the best fits with the MF result Eq. We see that 
all of the resulting distributions fit well with the MF result with \x ~ /-* m f- Figure [l] (b) 
is a case of regular networks with z = 0.001. We find that the results do not agree 
with the MF result, but instead fit well with log- normal distributions, which are shown 
by the solid lines. Note that the best fit log-normal distributions have time- dependent 
mean value mj n (t) and time-dependent variance cr? (t). On the other hand, if we rewire 
links with high probability the MF type distributions reappear. Figure [I] (c) is a case of 
regular networks with z = 0.001 and p = 0.5. The solid lines are the fitting by MF type 
solutions with \i = 1.7 ~ 1.8, which is slightly smaller than /i m f. Hence we observe that 
the change in wealth distribution occurs at the intermediate value of p. Figure [I] (d) is 
a case of small-world networks with z = 0.001 and p = 0.05. We have found that the 
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distributions can not be fitted well by either the log-normal functions or the MF type 
functions. Interestingly, however, we find that a combination of log-normal and power- 
law distributions fits well with the results. One such example at the 2 x 10 4 -th time step 
is shown with a log-normal function for middle and low wealth ranges (solid line) and 
a power- law function for the high wealth range (dashed line). Log- normal distributions 
with a power law tail like this are frequently observed in real world economic phenomena 
I- 

It is natural to question why distribution patterns change. To answer this question, 
we directly consider the development of the wealth. Examples of numerical simulations 
are shown in Figs. [2] (a)-(d). In these figures the horizontal axis shows the agents' number 
and the vertical axis shows time steps. Though we have simulated for N = 10 4 agents, 
wealth developments for 10 3 agents are shown. We have calculated up to 10 4 time steps. 
If any agent ranks in the top 10%, we mark it with a black dot, and if any agent ranks in 
the bottom 10% we mark it with a red dot. 

Figure ||] (a) is a case of regular networks with z — 0.1. We can see that the black and 
red dots are uniformly distributed without notable clustering. Since each agent has 10 3 
neighbors, each agent transacts with agents from various ranks, covering a wide range. 
Figure (b) is a case of regular networks with z = 0.001. Contrary to the case of z — 0.1, 
the black and red dots cluster heavily. This "wealth clustering" phenomenon is supported 
by the fact that each agent has only 10 neighbors and rich agents transact mainly with 
rich agents and poor agents transact mainly with poor agents. Hence Wi(t) and Wi(t) 
have strong correlations. Actually, for the case of regular networks with z = 0.001, we 
can obtain (Wi(t)wi(t)) = 0.92 ~ 0.95. Here w^t) = Wi(t)/W(t). On the other hand 
we can obtain (Wi(t)wi(t)) = 0.02 ~ 0.09 for the case of regular networks with z = 0.1. 
Figure |2] (c) is a case of small-world networks with z = 0.001 and p = 0.5. This case 
is almost similar to Fig. |^ (a). Figure ^| (d) is a case of small- world networks with 
z = 0.001 and p = 0.05. Although wealth clustering is observed, the size of each cluster 
is smaller than the case of regular networks with z = 0.001. From these results we find 
that a decrease in the number of links causes wealth clustering and the rewiring of links 
disperses it. 

If Wi(t) and Wi(t) have strong correlations, Eq. (^) becomes a pure multiplicative 
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Figure 2: Wealth developments for regular networks with z — 0.1 are shown in (a) and 
with z = 0.001 in (b). The case of small-world networks with z = 0.001 and p = 0.05 is 
shown in (c), and p = 0.5 in (d). 

stochastic process, and the wealth obeys log-normal distributions with mean vt and vari- 
ance Dt, where v = (\og{a(t) + a 2 }} and D = ([\og{a(t) + a 2 }} 2 } - (log{a(t) + a 2 }} 2 . 
On the other hand if Wi(i) and Wi(t) have no correlations, Eq. (^) is regarded as a 
multiplicative stochastic process with additive noise ||. It is known that this process 
induces power law distributions in large Wi(t) ranges, if two conditions are satisfied: One 
condition is that the stochastic variable and the additive noise JWi{t) must behave 
independently. The other condition is (log{eij(t) + a 2 — J}) < 0. Although the first con- 
dition does not matter in our model, the parameter ranges of J and a 2 are restricted by 
the second condition. If Wi(t) and Wi(t) have no correlations in the high wealth ranges 

6 



(a) (b) 




1 2 3,4 5 1 2 3 , 4 5 

Time Steps (1 ) Time Steps (10 ) 



Figure 3: Change of Gini coefficient in the case of regular networks (a), and that of 
small- world networks (b). 

and have strong correlations in the middle and low wealth ranges, the wealth will obey 
log-normal distributions with a power law tail such as those observed in many real world 
economic phenomena ||. 

We will now study the inequality of wealth distribution by using the Gini coefficient 
G to quantify it. Take every possible pair of income recipients (wealth possessors) and 
calculate the average of the absolute differences of the two incomes (wealth) for each pair. 
G is then half the ratio of the average to the mean value of all the incomes (wealth). 
Obviously G = when everyone has the same, and G = 1 when the entire amount is 
concentrated in a single person. With our parameter set, the Gini coefficient of the MF 
case is G m £ = 0.5. 

The temporal changes of G are shown in Figs. || (a) and (b). In these figures, the 
horizontal axis shows time steps and the vertical axis shows the value of G. We have 
numerically simulated up to 5 x 10 4 time steps. The temporal changes of G in the case 
of regular networks are shown in Fig. ^| (a). In this figure we observe that G gradually 
increases with time steps in the case of z = 0.01 and 0.02, which is in accordance with the 
log- normal distributions with mj n (t) and &f n {t)- We also observe that almost completely 
inequal distributions appear in the case of z = 0.001. We find that the decrease of z makes 
distributions inequal. The temporal changes of G in the case of small-world networks with 
z = 0.001 are shown in Fig. |3| (b). From this figure we find that the rewiring of links in 
small-world networks makes distributions equal. 
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Figure 4: Patterns of wealth distribution in the case of N = 10 4 . 

The results of this study are summarized in two parts. First, we have carried out a 
simulation for a wide range of values of p and z, which is summarized by the phase diagram 
shown in Fig. [|. This diagram has two extreme distributions, the log-normal distributions 
with mi n (t) and and the MF type distributions represented by Eq. @. These two 

ranges are interpolated by a decrease in the number of links or the rewiring of the links. 
In the intermediate region, the wealth distribution may be fitted well with log-normal 
distributions with a power law tail, which are observed in real world economic phenomena 
|1. Lastly, from the study of the Gini coefficient we have found that a decrease in the 
number of links makes wealth distribution inequal and the rewiring of links in small- world 
networks makes it equal. Hence society will be protected against exceeding inequalities 
in wealth by small-world effects in economic networks. 

We note that it would be very interesting to carry out an analysis of models similar 
to the ones studied here, but built on other kinds of networks. For example, the network 
can be dynamic, allowing rewiring to be either independent from or dependent on wealth. 
Other novel networks, like scale-independent networks, would also be worth studying. 
We believe that the work presented here forms a basis for the study of general stochastic 
processes on complex networks. 
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